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ABSTRACT 

We report monitoring observations of the T Tauri star EX Lupi during its outburst in 2008 in the 
CO fundamental band at 4.6-5.0 /im. The observations were carried out at the VLT and the Subaru 
Telescope at six epochs from April to August 2008, covering the plateau of the outburst and the fading- 
phase to a quiescent state. The line flux of CO emission declines with the visual brightness of the 
star and the continuum flux at 5 fim, but composed of two subcomponents that decay with different 
rates. The narrow line emission (50 km s _1 in FWHM) is near the systemic velocity of EX Lupi. 
These emission lines appear exclusively in v=l-0. The line widths translate to a characteristic orbiting 
radius of 0.4 AU. The broad line component (FWZI ~ 150 km s" 1 ) is highly excited upto v<6. The 
line flux of the component decreases faster than the narrow line emission. Simple modeling of the 
line profiles implies that the broad-line emitting gas is orbiting around the star at 0.04-0.4 AU. The 
excitation state, the decay speed of the line flux, and the line profile, indicate that the broad-line 
emission component is physically distinct from the narrow-line emission component, and more tightly 
related to the outburst event. 

Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: activity: 
other — stars: pre-main sequence — stars: formation — stars: individual (EX 
Lupi) 



1. INTRODUCTION 

EX Lupi, the prototype of the class of the eruptive 
pre-main sequence stars called EXors Q underwent its 
largest outburst recorded in the past 50 years in early 
200 8 (history of the ob servations of EX Lupi and EX- 
ors: IHerbid 120071 12008D . The observational definition of 
EXor variables relies strongly on the prototype EX Lupi, 
which spectroscopically l ooks like a classical T Tauri star 
in the quiescent phase (jHerbiel 119501 |2001[ ) , but flares 

mgoto@mpia.de 

1 Based on data collected in the course of CRIRES Large pro- 
gram [179.C-0151] and DDT program [281.C-4031(D,E,F)] at the 
VLT on Cerro Paranal (Chile), which is operated by the Euro- 
pean Southern Observatory (ESO). 

2 Based on data collected at Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan. 

3 Max-Planck-Institut fur Astronomie, Konigstuhl 17, D- 
69117 Heidelberg, Germany. 

4 Konkoly Observatory, Konkoly Thege Miklos 15-17, H-1121 
Budapest, Hungary. 

5 ESO, Karl-Schwarzschild-Strafie 2, D-85748 Garching bei 
Miinchen, Germany. 

6 Max-Planck-Institut fur extraterrestrische Physik, Giessen- 
bachstraBc D-85748 Garching bei Miinchen, Germany. 

7 ISDC Data Centre for Astrophysics & Geneva Observatory, 
University of Geneva, chemin d'Ecogia 16, 1290 Versoix, Switzer- 
land. 

8 California Institute of Technology Division of Geological 
& Planetary Sciences, MC 170-25 1200 E. California Blvd., 
Pasadena, CA 91125, USA. 

9 Johns Hopkins University, Department of Physics and As- 
tronomy 3701 San Martin drive Baltimore, MD 21218 USA. 

10 Leiden Observatory, Leiden University, P.O. Box 9513 NL- 
2300 RA Leiden, The Netherlands. 

11 Subaru Telescope, 650 North A'ohoku Place, Hilo. HI 
96720, USA. 

12 The gr oup of l ow-m ass stars with large variability is coined 
as EXors by Hcrbig (1989) for its nature of the outburst similar to 
FU Ori variables that are collectively called FUors. 



by 1-3 magnitudes in the visible roughly once a decade , 
with each outbreak lasting about a year (jHerbid 12007). 
The spectroscopic behavior in the outburst is compli- 
cated, an d sometimes related to an expanding stellar 
shell (e.g.|Herbig|[l989). It is increasingly accepted, how- 
ever, that the circumstellar disk plays the primary role in 
the outbursts, in which mass accreted from the envelope 
is stalled somewhere in the disk, until it reaches a critical 
surface density and triggers a disk instability. 

The idea that EXor outbursts arise inside the accre- 
tion disk originates in their observational similarity with 
FUors. The most compelling evidence that the FUor 
outburst has its origin in a disk is the split absorption 
lines in the visible to the near-infrared with wavelength 
separation of two absorpt ion peaks gradually decreas- 
ing w ith the wavelength (|Hartmann. Hinkle. fc Calvetl 
2004). This is most naturally accounted for when 
the absorption lines at long wavelengths arise from the 
cool outer part o f the disk , whe re the orbital veloc- 
ity is slower (e.g. iZhu et all 120091 ). FUors and EXors 
are historically considered as different classes of erup- 
tive pre-main sequence stars. However, both classes 
may represent a similar eruptive event, and differ only 
quantitatively in terms of flare amplitude, outburst du- 
ration, and outburst frequency. Recent outbursts of 
V1647 Ori (the central so urce of "McNeal's Nebula" , e.g. 
lAspin. Beck, fc Reipur th 2008) are recurrent like EX- 
ors, yet the amplitude of the outburst is more compat- 
ible with FUors, which fills the gap between these two 
classes. The mass accreted in a single outburst and the 
frequency of the outbursts indicates that an EXor/FUor 
star acquires a substantial fraction of its mass solely dur- 
ing the outbursts. It lends circumstantial support to 
the theory that EXors/FUors are not a group of spe- 
cial variable stars, but a state of the accretion phase 
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Fig. 1. — Light curve in the visible (crosses; from AAVSO http://www.aavso.org/ ) with the ordinate to the left, and in M band 
(diamonds for CRIRES, and a square for IRCS) to the right. The axis of M band magnitude is arbitrarily scaled to that of V band. The M 
band photometry is performed with the spectroscopic data, with the continuum level of the spectroscopic standard star as the photometric 
reference. The uncertainty of M band photometry is ±0.5 mag (the error bar shown in the upper right corner). 



that stars commonly experienc e during their evolution 
(e.g. lHartmann fc Kenvon[|1996D . However, there are in- 
deed qualitative differences between FUors and normal 
T Tauri stars, such as the clustering nature, as FUors 
tend to be in isolation, while other T Tauri stars are 
in clusters dHerbig. Petrov. fc Duemmlerl 120031 ). More- 
over, the number of known FUors in the Orion Neb- 
ula Cluster falls short if all class I source s go through 
the F Uor phase with a finite duty cycle (jFedele et al.l 
I2007T ) . although the number may increase with the in- 
clusion of the "FUor-like" stars b eing identified recently 
(|Greene, Aspin. fc Reipurthll2008l) . The current focus of 
EXor studies is to understand (1) what is the mecha- 
nism of the disk instability, (2) what is the relationship 
between EXors and FUors, and (3) what is the role of 
EXor/FUor events in the context of star formation. 

In this report, we will concentrate on the cause of 
the disk instability and the physical scale of the out- 
burst, as it often gives a clue to the trigger of the in- 
stability. The most favored model to date, thermal in- 
stability in a disk, is induced by ionization, and sub- 
sequent opacity increase due to negative hy drogen ions 
(|Bell. Lin fc Rudenlli"99l IClarke et all [19891) . The ther- 
mal instability is therefore only operational within a lim- 
ited disk radius, where atomic hydrogen is subject to 
thermal ionization. If the physical size of the outburst 
is larger, magnetohydrodynamic (MHD) instability may 
be favored, as it predicts a ma ssive accretion in the m ag- 
netic dead zone at 1-10 AU ([Armitage et al.l 120011) . If 
the outburst involves the whole disk, the scenario of an 
unseen co mpanions leading to globa l disk instability may 
also work (|Bonnell fc Bastienlll992T) . 

The vibrational transitions of CO are a unique probe 
of the inner regions of disks (0.01-1 AU) due to their 
high critical density (n cr > 10 10 cm~ 3 ) and high excita- 
tion energy (AP v= i_ > 3000 K). Measurement of the 
location of the emitting gas is relatively straightforward 
under the assumption that th e gas is in Keplerian ro- 
tation (e.g. INaiita et al.1ll996D . We utilized multi-epoch 
observations of CO vibrational transitions in EX Lupi to 
constrain where in the disk the hot gas is located dur- 
ing the outburst, and how its physical property devel- 
ops with time. The present study is part of the first 
multi- wavelength campaign of EX Lupi from the optical 



(|Sicilia-Aguilar et al.l 120101 ) . near-infrared ([Kospal et al.l 

20101) to mid- i nfrare d wavelengths ((Abraham et al.|[2009t 
Juhasz et al.1 [2010) during the 2008 eruption. We 
will use d = 155 pc as the distance to EX Lupi 
from ILombardi et al.l (I2008D , and the stellar parameters 
fromlGras- Velazqu ez fc Ravi (|2005l ) following lSipos et"al"1 
((20090 (Table [[]). 

2. OBSERVATIONS 

The spectroscopic observations of the CO fundamental 
band were performed at the VLT and the Subaru Tele- 
scope over a total of 6 epochs (Figure [1} . Half of the 
data were obtained in the director's discretionary time 
(DDT) on 14, 17, and 21 August 2008 at the VLT. Two 
wavelength intervals of the CO fundamental band (4.66- 
4.77 /Ltm and 4.86-4.99 /im with gaps betw een the four 
detec tor chips) were observed with CRIRES (Ka un et al.l 
I2004D to cover the low [P(l)-P(ll)] and the high [P(21)- 
P(32)] rotational transitions of v=l-0 in order to pro- 
vide good constraints on the rotational excitation tem- 
perature of the hot gas. The slit width was 0'.'2, corre- 
sponding to a nominal spectral resolution of i?=100,000 
(Aw=3 km s _1 ). The spectra were recorded by nod- 
ding the telescope along the slit to remove the sky back- 
ground. The slit was oriented along P. A. = 0°, 180°, 
90° and 270° by rotating the image rotator for the possi- 
ble use of the data for spectroastrometry. Unfortunately, 
the achieved signal to noise ratio was not sufficient to 
extract any meaningful astrometric signal. The data 
obtained with different P.A. were simply summed and 
reduced as classical spectroscopic data. The total inte- 
gration time was 12 and 16 minute for eac h grating set- 
ting. The MACAO adaptive optics system ([Bonnet et al.1 
120041 ) was used with EX Lupi itself as the guide star. 
Spectroscopic standard stars were observed at close air- 
mass with a same instrument setup. 

EX Lupi was also observed by CRIRES at the VLT in 
the framework of the ESO large program 179.C-0151 on 
April 28 and August 7 in 2008. The observations covered 
the CO spectrum at 4.65-4.77 jum with additional cov- 
erage at 4.81-4.90 /jm in the August observations. The 
slit width and the AO settings were the same as during 
the DDT observations. The total exposure times were 
about 20 min per grating setting in both of the nights. 
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TABLE 1 
Stellar Parameters for EX Lupi 

USED IN THIS PAPER. 

Reference 



M* ... 


0.6 M© 1. 


R* ... 


1.6 Rq 1. 


L, ... 


0.5 L© 1. 


d a . . . 


155 pc 2. 


V LSR b ■ • ■ 


f»0 km s" 1 3. 


Note. 


- Reference: 



1. 



iGras- Velazquez k Ravi (12005), 



2. ILombardi et all (120080. 

3. iHe rbig, Pe trov. k D ucmmlcr 
([200311. 

a Distance to EX Lupi. 

b Systemic velocity as measured 

with optical spectroscopy. 



IRCS (|Tokunaea et al.lll998l ) at the Subaru Telescope 
was used to obtain CO fundamental band spectra of 
EX Lupi on 21 July 2008. The angle of the echelle grating 
was set so that the spectral interval of 4.65-4.74 /jm was 
covered by the detector. The 0'.'15 slit was used resulting 
in a spectral resolution of i?=20,000 (Av — 15km s _1 ). 
The spectra were recorded by nodding the tip-tilt mir- 
ror inside the AO system 3'.'4 along the slit to remove the 
sky background. The observation was performed through 
a large airmass (secz ~ 2) with 5 times lower spectral 
resolution than CRIRES. Thus, the data quality is sig- 
nificantly lower. The total integrati on time was 14 min- 
utes. The adaptive optics system (|Havano et al.l I2008T) 
was also used here with EX Lupi as the guide star. The 
observations of all 6 epochs are summarized in Table [3] 

One-dimensional spectra were extracted from the 
CRIRES dataset using the crires_spec_jitter recip4___ on 
the ESO gasgano platform^. A custom-written IDL 
code was used to divide the spectra of EX Lupi by 
standard star spectra to eliminate telluric absorption 
lines. In the cases where the precipitable water vapor 
varied during the night, the water absorption was re- 
moved from both the object and the standard spectra be- 
forehand by dividing them by th e ATRAN atmospheric 
transmission model ([Lord! |1992|) . The IDL code fur- 
ther corrected wavelength offsets between EX Lupi and 
the spectroscopic standard spectra, slight differences in 
spectral resolution, and mismatches of the airmass. The 
wavelength calibration was performed by maximizing the 
cross-correlation between the observed spectra and the 
atmospheric transmission curve. The wavelength cali- 
bration error depends on the density of the telluric lines, 
and is less than a few pixels (~5 km s _1 ). 

The IRCS data were reduced in the same way, except 
that the IRAlP^l aperture extraction package was used 
to extract the one-dimensional spectra. 

EX Lupi is an M0 star in its qui escent phase wit h 
CO bandhead absorption at 2.3 (im ([Sipos et all 120091. 

13 CRIRES Pipeline User Manual VLT-MAN-ESO-19500-4406. 

14 http:/ /www. eso.org/sci/data-processing/software/gasgano/. 

15 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Nevertheless, the photospheric correction for the CO ab- 
sorption lines was not performed at 4.7 /xm, as the pho- 
tospheric temperature during the outburst is unknown. 
The contribution of the stellar photo sphere at 5 <nm i s 
less than 20% in the quiescent phase dSipos et alJl2 0091. 
and e xpected to be less during the outburst (Juha sz et all 
[20101) . 

Since there is no simultaneous photometry available at 
5 /im, the absolute flux calibration was performed against 
the spectroscopic standard stars HR 6084, HR 5987, and 
HR 6508. HR 6084 is a Bl III star with an M band 
magnitude of 2.42 ma£E HR 5987 and HR 6508 are 
B2.5 Vn and B2 IV stars with K magnitudes of 4.70 and 
3.18 mag, respectively, but without accurate M band 
brightness documented. We applied a color correction 
of K — M = — 0-05 m ag taken from K — L of early B 
stars iTokuna ga (1999) with the original from lKoornneeJ 
(1983)], that gives M = 4.75 and 3.23 mag, respectively. 
The precise correction is unimportant, since the overall 
photometric accuracy is only ±0.5 mag. The main source 
of photometric error is the variable slit transmission at 
the time of the observation up to a factor of three even in 
a pair of consecutive frames. The variable atmospheric 
transmission during the nights adds another 30% to the 
photometric uncertainty. EX Lupi becomes fainter from 
3.9 to 5.6 mag in the M band over 4 months (Table [J). 
The photometric value on 28 April is surprisingly bright, 
as Spitzer spectroscopy in A pril 2008 indic ates 1.5 Jy at 
5 /xm, or M «5.0 mag ([Juhasz et alJl2010h . Since then, 
the star became fainter at optical wavelengths^, albeit 
with significant variability on timescales of days to weeks. 
The data counts in the raw frames are at least consistent 
with those of M=3.9-5.6 mag stars, as simulated by the 
exposure time calculator of CRIREgP^l. 

3. ANALYSIS 

3.1. Two Emission Components 

All spectra from the 6 epochs are shown in Fig. [2] af- 
ter having been normalized to the continuum and cor- 
rected to the local standard of the rest. Two spec- 
tral components are immediately noticeable in the CO 
spectra, one broad [full width at zero intensity (FWZI) 
~ 150 km s _1 ] and one narrow [full width at half maxi- 
mum (FWHM)~ 50 km s _1 ] overlaid on each other. The 
broad line emission is characterized by the hot vibra- 
tional transitions up to v=6 confirmed. This is in agree- 
ment with the first overtone band in emi ssion seen at 
least up to v=5-3 in the SINFONI spectra ([Kospal et all 
120101 ) observed on 25-31 July prior to the CRIRES DDT 
runs during the outburst. In order to illustrate the broad 
line profiles of the former, a few lines from v>2 levels, 
which are relatively isolated from others, such as P(3) 
v=2-l at 4.750 ^m and R(9) v=3-2 at 4.703 /Ltm, are 
combined and shown in the left panel of Fig. [3] The 
broad emission lines are most clearly seen in the April 
28 data, taken during the outburst, a few months after 
the optical maximum. 

The sharp emission lines on top of the broad com- 
ponent are close to the systemic velocity of EX Lupi 

16 http:/ /www .jach.hawaii.edu/UKIRT/astronomy/calib/phot_ 
cal /bright_st ds.html 

17 AAVSO http://www.aavso.org/ 

18 http:/ /www. eso.org/observing/etc/. 
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TABLE 2 
Journal of observation. 



Epoch 



Program 



.ID 



V 
[mag] 1 



M Spectrograph Telescope Velocity 
b Resolution 



[mag- 



Wavelength Integration Standard Spectral 

Coverage Time [min] c Star Type 



28 Apr 2008 179.C-0151(D) 2454585 9.7 3.88 CRIRES VLT 3 km s" 

21 Jul 2008 DDT 2454669 10.0 4.50 IRCS Subaru 15 km s" 

7 Aug 2008 179.C-0151(D) 2454686 10.8 4.18 CRIRES VLT 3 km s~ 

14 Aug 2008 281. C-4031(D) 2454693 11.1 5.27 CRIRES VLT 3 km s" 

17 Aug 2008 281.C-4031(E) 2454696 12.3 5.53 CRIRES VLT 3 km s" 

21 Aug 2008 281.C-4031(F) 2454700 13.3 5.58 CRIRES VLT 3 km s~ 



1 4.65-4.77 jim 

1 4.65-4.74 ura 

1 4.65-4.77/im, 4.80-4.90Atm 
1 4.66-4.77/im d 4.86-4.99jtm d 
1 4.66-4.77/im d 4.86-4.99Atm d 
1 4.66-4.77/im d 4.86-4.99jtm d 



20 

11 

20 
12, 16 
12, 16 
12, 16 



HR 6084 
HR 5987 
HR 6084 
HR 6508 
HR 6508 
HR 6508 



Bl III 
B2.5Vn 
Bl III 
B2IV 
B2IV 
B2IV 



a From AAVSO [http: / / www. aavso.org/ 

k Photometry is done with respect to the spectroscopic standard stars. The accuracy is ± 0.5 mag. 
c For one grating setting. 

d With gaps between the four detector chips. 



Ksr ~ km s" 1 ; IHerbig. Petrov. fe Duemmlerl 120031) . 
The intensity of the narrow lines is 10-20% above the 
continuum level throughout the observations. The line 
profile of the emission lines does not change significantly 
during and at the end of the outburst (Fig. [3l right). 
The line widths, markedly smaller than the broad line 
component, imply that the emission arises from regions 
of the disk outside where the broad and highly-excited 
emission lines comes from. 

The spectral profile, composed of two distinct kine- 
matic components, is somewhat similar to Fe II at 
5018 A observed with the HIRES spectrograph (R = 
48, 000) at Keck in 1998 when EX Lup i was at the max- 
imum in a previous outburst. iHerbia (|2007l ) attributed 
the broad line emission to the accretion funnels, as the 
Fe II emission had inverse P Cyg profile with broad ab- 
sorption longward at +300 km s _1 . The velocity widths 
of the broad and the narrow components are 140 km s _1 
and 23 km s _1 in Fe II 5018 A, similar to the CO lines 
at 4.7 /mi. We looked for similar systematic redward ab- 
sorption in the CO lines, and found that they are not 
present in either epoch of the observations. Moreover, as 
we discuss later, the line profile in the outburst is double- 
peaked, which implies that the emitting gas is orbiting 
around the star. We conclude that the origin of the broad 
CO emission is not the stellar photosphere or accretion 
columns, but the circumstellar disk. We will deconvolve 
the broad and the narrow emission components using a 
simple slab model in the next sections. 

3.2. Population Diagrams 

First, population diagrams were constructed to obtain 
an overview of the physical property of the hot emit- 
ting gas in the disk. Equivalent widths of the emission 
lines were measured integrating v=l-0 P-branch lines 
over ±120 km s _1 . They thus include v=l-0 emission 
of both broad and narrow components without correct- 
ing for the overlapping higher vibrational transitions. 
The column density at each rotational level is calcu- 
lated as huNjA v ■ ttR 2 = A-Kd 2 W u fo from the equiva- 
lent widths W v . fo is the continuum flux coming from 
the M-band photometry. The specific radius of the emit- 
ting area R s is set to 0.2 AU, so that it is roughly con- 
sistent with the line width of the broad line emission. 
The column density normalized by the degeneracy of the 
levels Nj/(2J + 1) is plotted against the energy levels 
Ej/k in Fig. [4j The observed population diagrams show 



turnovers at Ej/k — 3200-3400 K, which testifies the 
emitting gas is optically thick. The vertical offset of the 
diagrams over 6 epochs spanned roughly an order of mag- 
nitude. The population diagram was initially fit by a 
thermodynamical model with T x and iVco as free pa- 
rameters. However, in order to reproduce the sequences 
of Nj in the optical thick regime with an order of mag- 
nitude dispersion, Nqo has to be adjusted by several or- 
ders of magnitude up to 10 25 cm" 2 , which is apparently 
overstretching the capability of the model. We there- 
fore fixed the column density to the minimum value of 6 
epochs (AVjo = 2 x 10 20 cm~ 2 ) found in the first path 
of the fitting, and tried to match the absolute levels of 
Nj by changing the size of the emitting area. The col- 
umn density 2 x 10 20 cm~ 2 is close to the upper limit 
of the physically feasible value of the relatively small 
area of the disk {R s ~ 0.1 AU, H/R=0.1, [C/H] = 10~ 4 , 
tih = 10 12 cm -3 ). The fractional column density at v=l 
level was calculated with the vibrational temperature T v 
set equal to T x . The best-fit T x is constrained within 
1100-2400 K and tabulated in Table [3] with the corre- 
sponding emitting radii of the disk R s . The line profile 
is assumed to be a Gaussian function of o~ v = 5 km s _1 . 

3.3. Slab Model 

Second, we used a simple slab model to incorporate 
the geometry and the rotation of the disk to repro- 
duce the observed spectra. The slab model is a gas 
disk with no geometrical thickness, and meant to re- 
produce the underlying broad line emission only. The 
fundamental transition of CO is used as the sole opac- 
ity source at the wavelength concerned. No isotopomers 
are included, as this does not change the quality of the 
fit within the given signal-to-noise ratio of the spec- 
tra. Rotational levels up to J = 100, and vibrational 
levels up to v=8 are taken into account. The wave- 
lengths, en ergy levels, and E instein A v coefficients are 
taken from iGoorvitch I (|1994h . Ro- vibrational levels are 
populated according to the single-temperature Boltz- 
mann distribution. Thermal broadening at the rotational 
temperature is taken as the only source of the intrin- 
sic line width. The input parameters of thermodynam- 
ics, the total column density Nqoi and the rotational 
and vibrational temperatures T x and T v , are kept con- 
stant without radial dependency. The resulting spec- 
trum is calculated as I v = -^S v {\ — e _T "), where S v 
is the line source function A2iri2/(-Bi2^i — #21^2)1 and 
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Fig. 2. — Spectra of EX Lupi at the CO fundamental band over six epochs from 28 April to 21 August 2008. All spectra were obtained 
by CRIRES at the VLT (R= 100,000) except on 21 July, when IRCS at Subaru Telescope was used with 5 times lower spectral resolution 
(i?=20,000). The spectra are normalized to the continuum level. Two emission components are immediately noticeable: one broad emission 
that fades away in August, and narrow line emission. Two model spectra for the outburst component are shown overlaid. The blue models 
are calculated with the excitation temperatures that fit the population diagrams in Fig. |4]and the constant column densities 2 X 10 20 cm -2 . 
The red models are calculated by fixing the excitation temperature at T x = 4500 K, which is the dissociation limit of the molecules. The 
two models supposedly represent the extreme cases meant to illustrate how the thermodynamical parameters are degenerated. The model 
parameters are listed in Table [3] The spectra are not covered continuously because of the gaps between the detector chips in CRIRES. 
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Fig. 3. — Left: composite line profiles of the vibrationally excited 
transitions (v>2). The gray lines are the overlay of CO P(3) v=2-l 
4.750 Atm, CO P(4) v=2-l 4.759 urn, R(9) v=3-2 4.703 um, and 
-R(21) v=4-3 4.679 fim. The black line is the averaged spectrum. 
The line emission is double-peaked. Right: composite line profiles 
of CO v=l-0 after the outburst models of T x = 4500 K in Figured 
are subtracted. The gray lines are overlay of v=l-0 lines within the 
spectral coverage, and the black lines are the averaged spectrum of 
them. A blue depression shows up at the velocity of —80 km s 
in August, and becomes deeper later. We attribute it to the disk 
wind started at the late phase of the outburst. 

t v = ■^ 2 (n 1 i3 1 2 — 712-821) 4>vl cosi is the opacity with </>„ 
being the dimensionless line profile, and i being the disk 
inclination angle. 

We tested a model with the opacity of the dust grains 
uniformly mixed with the gas. If we take the mass ab- 
sorption co efficient k(X = 5 um) =2000 cm 2 g _1 , as cal- 
culated by Osscnkopf & Hcnning (1994), the line emis- 
sion of CO will be entirely suppressed by the dust con- 
tinuum emission, unless the gas to dust mass ratio is 
larger than 10,000. We infer that the continuum emis- 
sion at 5 /xm comes from the region outside the broad line 
emission region. The continuum emission is therefore re- 
moved from the observed spectra, and only the remaining 
line emission part is produced by the slab model, without 
dust continuum opacity in the radiation. 

The disk is divided into annuli from the inner (R4) to 
the outer disk (R ) boundary in the way that logr,-+i — 
logrj is constant. The ro- vibrational emission spectrum 
is calculated in a small segment of each annulus with 
the given rotational and vibrational temperatures. The 
line spectra are then shifted by the radial velocity of 
the Keplerian rotation at the given azimuthal location, 
and added toge ther. The stellar mass is a ssumed to 
be M* =0.6M Q (jGras- Velazquez fc Ravi [2005). The line 



emission spectra at the different radii are combined with 
an arbitrary weighting function I(r) = I(Ri)(-^-)~ a , 
where a is fixed to 2.5 in the present modeling. The 
inclination angles between 40° and 50° give reasonably 
good fits to the spectra at all the epochs. We there- 
fore fixed it t o i =45° wit h out fu rther adjustment. The 
disk model of iSipos et al.1 ((2009), which reproduces the 
optical-infrared SED of EX Lupi in the quiescent phase 
favors face-on geometry wit h the most plausib le incli- 
nation angle being i =20°. ISipos et al.l (|2009l ) found, 
however, that changing the inclination angle from 0° to 
40° changes the output SED little; as the SED is rela- 
tively insensitive to the inclination angle, unless the disk 
is heavily optically thick. The inclination angle we adopt 
here is near the upper limit of this range. 

4. DISCUSSION 
4.1. Spectral Components 
4.1.1. Broad Line Emission 

The disk parameters of the broad line emission that we 
modeled in the previous section, i, Ri, Nqo, T x and T v 
are not independent, but interlocked in a non-trivial way. 
The line width of vsini ~ 150 km s -1 on Apr 28 testi- 
fies that the radius of the inner edge of the disk is not 
larger than 0.09 AU. On the other hand, the disk can- 
not be too small in order to account for the absolute line 
flux. In the cases when the line profile is not simple and 
the signal-to-noise ratio is not high, such as the present 
case, the photometric calibration is important because 
the absolute line flux effectively constrains the size of 
the emitting region and therefore the disk radius. If we 
take the photometry at 5 /<m derived based on the spec- 
troscopic data as face value, the single line flux ranges 
from 10 -14 to 10 -13 erg cm -2 . The maximum flux den- 
sity that we could retrieve from a disk at d pc away is 
/„ « 7rr 2 i3„(Tb) / 'd 2 , where B v is the Planck function. 
Given the blackbody temperature of Tb = 2000 K, the 
distance of d = 155 pc, and the line width of 100 km s _1 , 
we will get f v < 10~ 14 erg s -1 cm~ 2 , if we assume a disk 
radius of 0.01 AU. This is too small to account for the 
observed line flux, therefore the radius of the inner hot 
disk cannot be smaller. The lower limit for the disk ra- 
dius then calls for the inclination angle to be larger than 
i >30° to match the observed line width. 

We first applied T x and iVco to the models as deter- 
mined by the optically thick population diagrams, and 
optimized gas radius to fit the line profiles best. The 
vibrational temperatures T v are set equal to T x . The re- 
sults underestimate the absolute line flux, because the 
outer radius R that is consistent with the kinematics in 
the line profiles is larger than the specific emitting radius 
R s we used to reproduce the population diagrams. The 
model spectra are normalized by the square of R a /R s 
and shown in Fig. [2] in blue lines. The absolute line flux 
is reproduced without such an arbitrary scaling, by in- 
creasing the excitation temperature, allowing it to differ 
from T v . However, when the emitting gas is hotter than 
3000 K, the line ratios of different rotational transitions 
diminish in a short wavelength interval. It becomes less 
straightforward to determine T x uniquely. In the opti- 
cally thick regime, increasing Aco and T x imposes sim- 
ilar effects on the output spectra. In order to cope with 
the degeneracy, and avoid increasing the optical depth 
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TABLE 3 

Disk parameters from our modeling. 



Epoch 


K a 
[AU] 


AT' b 
CO 

[xlO^cm- 2 ] 


rp! c 

[K] 


i d 
[°] 


Ri e 
[AU] 


Ro^ 
[AU] 


a s 


N CO h 
[xl0 19 cm- 2 ] 


Tv 1 
[K] 


Tj 
[K] 




Fx 1 

[ xlO" 13 W m-Vm- 1 ] 


28 Apr 2008 


0.5 


20 


2100 


45 


0.05 


0.3 


2.5 


6 


3000 


4500 


0.7 


5.9 


21 Jul 2008 


0.3 


20 


2400 


45 


0.05 


0.3 


2.5 


1 


2500 


4500 


1.0 


3.4 


7 Aug 2008 


0.4 


20 


2200 


45 


0.055 


0.4 


2.5 


2 


2700 


4500 


0.8 


4.5 


14 Aug 2008 


0.4 


20 


1300 


45 


0.05 


0.3 


2.5 


1.5 


2000 


4500 


-0.6 


1.7 


17 Aug 2008 


0.5 


20 


1100 


45 


0.04 


0.3 


2.5 


4 


1800 


4500 


0.6 


1.3 


21 Aug 2008 


0.4 


20 


1200 


45 


0.04 


0.2 


2.5 


2 


2000 


4500 


-0.8 


1.2 



a Radius of the emitting area from the population diagram in Fig. ^ 
b Total column density of CO fixed to 2 X 20 cm -2 . 

Rotational excitation temperature from the population diagram in Fig. ^ 
c Inclination angle from face-on. 
^ Inner radius. 
e Outer radius. 

^ Weight function of the intensity at different radius I(r) — I (Ri)(-£-)~ a - 

s Total column density of CO. 

k Vibrational excitation temperature. 

1 Rotational excitation temperature, fixed to 4500 K. 

J Asymmetry factor of line profile. 

k Flux density at the continuum level. 

any further, which our simple model cannot handle, the 
rotational temperature was set to 4500 K, as the thermal 
disso ciation temperature of CO is 4500-4600 K (|Tatuml 
1966). T v is adjusted so that the model reproduce the 
absolute line flux under the condition of T x = 4500 K. 
The second model with T x =4500 K is shown in Fig. [2] in 
red lines. The two solutions (Table |3|) can be regarded as 
the two extreme cases of the spectral models, and provide 
a fair idea how the model is affected by the degeneracy 
of the parameters. The best-fit line profile of a disk in 
Keplerian rotation sets the inner and the outer radius at 
0.04-0.05 AU (=5-7 R*) and 0.2-0.4 AU, respectively, 
for all epochs with the fixed inclination angle of i =45°. 
These geometric parameters are commonly used in the 
two representative models discussed above, and less sub- 
jected to the degeneracy. The optical depth is t u =0 
1000 at the line center in the model with T x — T v , and 
r„ =0-100 for T x = 4500 K, assuming the line is only 
thermal broadened. 

Although the heavy overlap of the lines makes it 
difficult to discern the shape clearly, the line profile 
looks double-peaked (Fig. [5]). This is the first time 
that the line profile is found to be double peaked in 
EXors. It lends support to the underlying similarity 
of the two classes of variable stars, FUors and EX- 
ors, and the origin of EXor outbursts associated with 
the circumstellar disk. Hartmann and his school (e.g. 
Hartmann fe Kenvo n 1996; Hartmann, Hinkl eT^g Calvetl 
20041 : iZhu et aLll2009l) attribute the split absorption lines 
of FUors, unusually broad for a stellar photosphere, to 
the gas in Keplerian rotation in unstable disks. The in- 
terval of line splitting is smaller in the near-infrared than 
in the visible. This is naturally explained if the absorp- 
tion lines at longer wavelengths a rise from the co oler, 
more distant, sl ow-rotating region (IZhu et al.ll2009l) On 
the other hand. | Herbig. Petrov. fc Duemmlerl (|2003f) and 
iPetrov fc Herbid (|2008l ) found no correlation between the 
line widths and the excitation potentials of the visible ab- 
sorption lines of FU Ori, which should be present if the 
lines arise from the disk, and the disk has a tempera- 
ture gradient. They argued against the disk origin of the 



optical lines, and claimed that the split lines are better 
reproduced by a large cool spot at the polar region on a 
rapidly rotating star. Moreover, the absorption lines of 
FUors are not perfectly Ke plerian, but have a rectangular 
shape at the bottom (e.g. lHartmann. Hinkle. fe Calvetl 
2004) . It is hard to distinguish whether the unusual line 
shape represents the turbulence in the upper layer of the 
disk through which the absorption takes place, or a par- 
ticular morphology an d geometry of the cool spot on the 
surface of a star (e.g. IPetrov fc H crbig 2008). The line 
profile of EX Lupi is found to be double peaked, but in 
emission. The double peaked emission lines would not 
be compatible with the stellar spot hypothesis, unless 
the stellar surface is extended to at least a few times the 
stellar radius. 

The line profiles are asymmetric between the blue and 
the red peaks. Figure [3] shows the blow-up view of the 
spectra near R(6) v=2-l (4.6675 /im). There are mi- 
nor contributions from i?(14) v=3-2 and i?(23) v=4-3 at 
4.6665 /im and 4.6670 /im, respectively. The asymmetric 
line profile, however, cannot be fully explained by chang- 
ing the vibrational temperature, thus the relative contri- 
bution of the three lines. The line profile is fitted better 
if we artificially introduce intrinsic asymmetry in the Ke- 
plerian line profile by modulating the azimuthal intensity 
distribution sinusoidally \I' V {&) — I V {9){1 + (1 — a) -sin 6*), 
where a is the asymmetric factor] by down to a =0.6. 
Asymmet ries in infrared C O line emission have been 
found by iGoto et al.l (|2006l ) in HD 141569 A, where the 
v=2-l emission is three times brighter in the northern 
disk than in t he southern disk in the sp atially resolved 
spectroscopy. iPontoppidan et al . (2008) proposed that 
the transition disks around SR 21 and HD 135344B have 
azimuthal structures, as the spectroastrometric signal of 
CO v=l-0 is asymmetric, which is consistent with th e 
millimeter dust continuum imaging (| Brown "eTaI1l2009T) . 
The line asymmetry of EX Lupi varies with time. 

4.1.2. Narrow Line Emission 

The single peak profile of the narrow line emission with 
50 km s _1 line- width in FWHM defies pure Keplerian ro- 
tation modeling, if the line intensity decreases as a power 



Goto et al. 




Apr 28 2008 
2100 K 



3000 4000 5000 6000 




Jul 21 2008 
2400 K 



3000 4000 5000 6000 



Aug 7 2008 
2200 K 




3000 4000 5000 6000 



Aug 14 2008 
1300 K 




10' 



.10' 



10' 



10' 



.10' 



10' 



10' 



.10' 



10' 



10' 



.10' 



10' 



10' 



10' 



10' 



10' 



.10' 



10' 

3000 4000 5000 6000 
EjA [K] 

Fig. 4. — Population diagram of CO v=l-0. The line flux is 
calculated by integrating over ±120 km s" 1 of the emission lines. 
The emission is optically thick as is shown by the turnovers at 
Ej /k = 3000—3400 K. The thermodynamical models with assumed 
line profile (Gaussian profile with a v = 5 km s —1 ) are shown over- 
laid. 



law of the radius. Instead, the line width is simply taken 
as the orbital velocity of the gas, and is translated to the 
disk radius of 0.4 AU, using the same inclination angle 
as the broad line emission. FWHM was used instead of 
FWZI, as the latter is hard to measure unambiguously 
because of the overlap with the broad line components. 
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The disk radius above does not mean the inner trunca- 
tion, but a characteristic radius where the narrow line 
emission arise from. The quiescent line emission comes 
exclusively from the v=l-0 transition with no obvious 
higher vibrational transitions. Except the emergence of 
absorption at the blue shoulder implying a disk wind, 
no significant changes are found either in the equivalent 
width or in the line shape (Fig. [6] Fig. [3l right). Note 
that the small equivalent width on July 21 is of less signif- 
icance than other epochs, as the data suffer from the low 
spectral resolution and the limited signal-to-noise ratio, 
which makes the removal of the telluric CO absorption 
challenging. The systematic error brought in by the tel- 
luric correction is magnified by the subtraction of the 
outburst model. The broad line emission is less subject 
to such problems in the measurements of the equivalent 
widths. The equivalent width of the narrow line compo- 
nent looks constant in time. 

4.1.3. Disk Wind 

The model spectra of the broad line emission system- 
atically overestimate the flux of v=l-0 lines in the blue 
shoulder, and underestimate the red. After subtract- 
ing the spectral model, therefore, the v=l-0 lines sys- 
tematically show a blue depression between —40 and 
— 100 km s _1 with excess emission at the red shoul- 
der at +80 km s^ 1 , which is similar to a P Cyg profile 
overlaid on the sharp quiescent emission (Fig. [31 right). 
This blueshifted absorption is more prominent toward 
the later epoch. The absorption lines can be noted at 
f (21) to P(26) on August 21, implying that the absorb- 
ing cloud is warm (T x ^100 K). The absorption is not 
likely photospheric, although EX Lupi is an M-type star 
in its quiescence, as the amount of blueshift is too large. 
The emergence of blue shifted absorption is also observed 
in infrared CO spectra of V1647 Ori toward the end of 
an outburst (|Brittain et al.ll2007f ) . We infer that the blue 
absorption in the EX Lupi spectra also represents the 
disk wind, similarly to the case of V1647 Ori. 

4.2. Heating Mechanism 

The excitation mechanism of the CO vibrational band 
is somewhat ambiguous. The elevated accretion energy 
can be directly converted to the kinetic energy of the gas 
through the viscous heating. In this case, the broad line 
emission that we locate at 0.04-0.4 AU traces the region 
where the accretion rate is high. On the other hand, the 
connection between the CO emission and the outburst 
could be indirect. The radiation from the central star, 
either from the hot accretion funnels or the footprints of 
the magnetic field lines, is first raised by the high disk 
accretion. The gas in the disk is heated afterwards by 
the enhanced irradiation. The primary difference of the 
two mechanisms is the vertical structure of the tempera- 
ture that the heating imposes. In the case that the gas is 
viscously heated by the accretion, the disk is hottest near 
the mid-plane, and gradually cools down toward the sur- 
face. The emergent spectr um should be in absorption, 
as is the case for FUors (jHartmann. Hinkle. Sz Calvetl 
2004). When the irradiative heating is dominant, the 
temperature gradient is inverse wit h emission lines from 
the upper layer of the disk surface (|Calvet et al.lll991ab . 

The irradiative heating has a few advantages to ac- 
count for the CO line emission of EX Lupi. First, most 
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Fig. 5. — Line profile of the R(6) v=2-l transition at 4.6675 fim. 
The systemic velocity of EX Lu pi is close to km s — 1 
IjHerbig. Petrov.~fc Ducmmlcr 20(53). The model spectra of the 
outburst with T x = 4500 K are overlaid in gray lines. There are mi- 
nor contributions from the higher vibrational transitions of -R(14) 
v=3-2 and i?(23) v=4-3 as shown in the right bottom panel on 
21 August. The observed line profiles cannot be reproduced with- 
out introducing the asymmetry. The asymmetry flips three times 
in August, which is consistent with the orbital period at 0.04- 
0.06 AU. 

of the infrared CO emission lines observed today are 
thought to come fro m the hot surface of a disk, irra- 
diatively heated (e.g. INaiita et all 120031) . This is fully 
consistent with the line flux of the narrow line emission 
closely following the continuum level. A critical problem 
of the direct heating of gas by the accretion is the cool- 
ing time of the outburst. The duration of the EX Lupi 
outburst as is seen in the visible light curve, is roughly 8 
months, from January to August 2008. The viscous time 
scale on which the gas accretes on the star is given by 

-L = 3 7 ra(H/r) 2 • -L 

^vis ^orbit 

for standard a prescription of an accretion disk. The vis- 
cous time scale at 0.04 AU is 150 to 600 days, depending 
on H/r from 0.05 to 0.025 at the radius, even for the 
fastest case with a being unity. The observed time scale 
of the outburst is un comfortably short c ompared to the 
most optimistic case ()Juhasz et al.|[20Toh . 

On the other hand, some features of the EX Lupi out- 



burst are not fully explained by the irradiative heating 
either. The line flux of both broad and narrow com- 
ponents decline with the visual magnitude and the con- 
tinuum flux at 5 /jm (Fig. [51 upper panels) , but on the 
slightly different decay rates. The narrow line emission 
is diminished by the factor of 4 from April to August 
(Table 2]) , which is the same factor with the decline of 
5 /xm continuum flux, while the broad line emission is 
decreased faster by the factor of 10. 

In order to contrast the decay rates, the equivalent 
widths are calculated in the same manner with the line 
flux. While the equivalent width of the narrow line emis- 
sion maintains the level of 1.5 x 10~ 4 jam at the end of the 
outburst (Fig[6l bottom right), where the visible magni- 
tude of EX Lupi returned to V =12-13 mag, compara- 
ble to the level of the pre-outburst brightness in 2007 
(Fig. [TJ), the equivalent width of the broad line emission 
declines, keeping pace with the continuum flux at the vi- 
sual wavelength. The differential time scale of the two 
line components is not accounted for, if the broad and 
the narrow line emissions are a single disk component 
radiative ly heated in the identical manner. 

In addition, the excitation state and the kinemat- 
ics are distinct in two emission components, with no 
smooth transition between them. The characteristic ra- 
dius where the narrow line emission arises from is 0.4 AU, 
and is constant as seen in the stable line width. The 
radius is immediately outside of the broad line emit- 
ting region at 0.04-0.4 AU (Fig. [7]). The location of 
this transition zone is close to the inner-rim of the dus t 
disk found by the SED modeling bv lJuhasz et all (|2010j) . 
With the vibrational excitation state upto v=6, and the 
faster decay of the line flux than the radiatively heated 
narrow line component, the broad line emission appar- 
ently traces the outburst event more directly. 

The heating mechanism of similar outburst variables, 
FUors, is most likely the elevated accretion near the disk- 
midplane, as the br oad and double-peak CO lines in ab - 
sorption testify (e.g. lHartmann. Hinkle. fc Cal vct 2004). 
To the contrary, the CO vibrational lines of EX Lupi 
are o bserved in emission in both overtone (Kosp al et al.l 
120101 ) and fundamental bands. The critical difference 
that might explain the emergent spectra of the two 
classes of the outburst variables is the optical depth of 
the disk. The broad line emitting region of EX Lupi 
is likely fre e of dust as i s disc ussed above. The SED 
analysis by Uuhasz et a l. (2010) also excludes the opti- 
cally thick continuum emission either by gas or by dust 
grains within 0.3 AU during the outburst. Moreover, the 
accretion rate during the outburst is a few orders magni- 
tude smaller in EX Lupi (~ 1O~ 7 M yr _1 ) than FUors 
(~ 1O~ 4 M0 yr -1 ), implying the lower disk surface den- 
sity. The disk close to the midplane that is viscously 
heated by the accretion might be more exposed to the 
surface in the case of EX Lupi, which possibly makes the 
CO vibrational band observed in emission. 

4.3. Origin of Outburst 

Let us extend our speculation to the trigger of the out- 
burst, taking the physical scale of it being traced by the 
broad line emission. The physical size of the outburst 
on the scale of 0.2-0.4 AU argues a gainst mechanisms 
that p redict a global disk instability. iBonnell fc Bastienl 
(1992) proposed that a binary system can tidally disrupt 
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TABLE 4 
Line flux and equivalent width. 



Epoch JD F\(V) Line flux Equivalent width 

[xlCT 12 Wm- 2 /im _1 ] [xlCT 17 Wm" 2 ] [xl0" 4 ^m] 

v>2 a v=l-0 b v>2 v=l-0 



28 Apr 2008 2454585 4.5 

21 Jul 2008 2454669 3.4 

7 Aug 2008 2454686 1.6 

14 Aug 2008 2454693 1.2 

17 Aug 2008 2454696 0.4 

21 Aug 2008 2454700 0.2 



9.7± 1.8 8.3± 0.8 1.6± 0.3 1.4± 0.1 

4.9± 2.5 2.7± 0.6 1.4± 0.7 0.8± 0.2 

5.5± 1.0 7.5± 0.5 1.2± 0.2 1.7± 0.1 

1.7± 0.5 2.3± 0.2 1.0± 0.3 1.4± 0.1 

1.0± 0.5 2.1± 0.2 0.8± 0.4 1.6± 0.2 

0.9± 0.4 2.2± 0.2 0.7± 0.3 1.7± 0.1 



a Integrated over ±120 km s 1 . 

^ Integrated over ±40 km s — 1 , after the spectral model with T x — 4500 K is subtracted. 




5 May Jun Jul Aug 

Flux Density at V [x10"' 2 W/m 2 /(m] [Dote] 



Fig. 6. — Top: (Left) Line flux of the emission lines from vibra- 
tionally excited states (v>2, open circles), and the lowest transi- 
tion v=l-0 (triangles) after subtraction of the outburst model of 
T x = 4500 K shown as functions of the continuum flux of EX Lupi 
at the visible wavelength. The line intensity is integrated between 
— 120 and +120 km s _1 for the broad line component. The com- 
posite line profiles used in the calculations are shown in Fig. [3] 
(Right) Same with the left panel, but shown with the date in ab- 
scissa. Bottom: Same with the top panels but for the equivalent 
widths. The line flux of the v>2 and v=l-0 lines decrease with the 
time and the visual brightness of the star, but with different speed. 
The difference is clearly seen in the plots of the equivalent widths. 
The equivalent widths of the broader, highly excited lines decays 
quickly with the optical brightness of the star, while the equivalent 
widths of the v=l-0 narrow line emission stays constant during the 
outburst. 

the disks of each other in a close encounter, and tem- 
porarily increase the accretion rate up to 10 M yr _1 . 
The tidal disruption by a companion is attractive because 
it naturally explains the repetitive nature of EXor out- 
bursts. The binary fraction among low-mass pre-mai n 
sequence stars is indeed high (e.g. iGhez et all [1993). 
FU Ori, the prot otype FUor, has a close companion 
(jWang et al.l [200 4V and there are even FUor- like bina- 
ries where both of the binary components are FUor-like, 
which is statistically unli kely if two are randomly paired 
(|Reipurth fc Aspinl I2004D . The gravitational disruption 
by a binary companion may work for the present case 
as a trigger of the disk instability combined with other 



mechanisms, but most likely does not explain the small 
radius of the high mass-accretion re gion by itself. No 
visible (IGhez et al.l 1 19971 iBailevi I1998D or spectroscopic 
companion (|Herbial2007D to EX Lupi has been identified 

a s of today. 

IVorobvov fc Basul ([20060 argued that planetary cores 
- which possibly have already formed in a gravitation- 
ally unstable disk - can trigger an FUor outburst when 
they migrate inward and eventually fall into the central 
star. Protoplanetary clumps form at > 10-50 AU away 
from the star. The migration from the outer disk over 
this distance does not explain the discrete boundary at 
0.2-0.4 AU either. Their simulation ends at the inner 
radius of 0.5 AU. Further investigation how the infalling 
a protoplanetary clump behaves in the final 0.5 AU is 
awaited. 

When mass inflow from the outer disk is faster than 
what the inner disk can transport further in, the accre- 
tion slows down, and piles up material at the boundary 
between the inner and the outer disk. The elevated sur- 
face density triggers a local disk instability, and the ac- 
cretion mode switches to a higher state. The mass accre- 
tion could slow down at a certain radius simply because 
the disk viscosity is too low, or the viscosity is locally too 
low because of the decoupling of the disk and the mag- 
netic field, or if the disk has a gap, or if it is magnetically 
truncated. 

The disk viscosity is provided by the coupling of the 
magnetic field and slightly ionized medium in the MHD 
instability model. Deep inside the disk where the stellar 
radiation no longer penetrates, there is a low ionization 
region called the dead zone, where the viscosity is close 
to zero. The incoming mass piles up there, until the 
surface density becomes high enough to heat the disk to 
«1000 K to restore the coupling with the magnetic field. 

A dead zone starts at 0.1 AU, and extends beyon d 
1 AU in a disk around a solar mass star (|Gammi e 199G). 
The disk instability is triggered at the c oldest region of 
the d e ad zone near the outer boundary ( Armitag e et al.l 
120011) . iZhu et aLl (|2007) constructed a radiation transfer 
model of the disk around FU Ori in outburst, and cal- 
culated the size of the active region to be around 1 AU, 
which naturally fits with MHD disk instability. The ob- 
served transition region of EX Lupi is a few times smaller 
than the typical size of a dead zone. This may not be 
surprising, as the luminosity of EX Lupi in quiescence 
is small as well [L* < 0.5Lq in Gras- Velazquez & Ray 
( 2005), L* = Q.lbL^ in Sipos et al. (2009)]. 

iWunsch et al.1 ()2005l ) discovered another type of insta- 
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Fig. 7. — Schematic representation of the broad and the narrow line emitting regions discussed in the text. 



bility near the inner edge of the dead zone in their nu- 
merical model. When a small perturbation is applied to 
the height of the dead zone, the mass inflow locally in- 
creases while the outgoing mass decreases. The height 
of the perturbed region grows, until it splits up from the 
dead zone and triggers a ring instability, or the inner edge 
of the dead zone starts oscillating radially (jWunsch et al.l 
2006). This instability works at the correct physical scale 
of 0.1-0.2 AU with small accretion rates 10 _9 M Q yr" 1 
which are appropriate for EX Lupi. However, the pre- 
dicted timescale of the outburst is 100 yrs or longer. This 
is close to the observed timescales of FUor outburst, but 
is too long for EX Lupi (~1 yr). 

When a star is strongly magnetized, the gas infall stops 
at the radius where the Keplerian angular velocity is 
equal to the stellar rotation, as the orbital motion of 
the gas is tied to the star. The matter piles up at the co- 
rotation radius, until the surface density becomes so high 
that the therm al pressure of the gas over comes the mag- 
netic pressure (|D'Angelo fc Spruitl 12010). This mecha- 
nism requires the presence of a strong magnetic field. 
Among young eruptive stars, FU Ori is th e first where 
a ma gnetic field was discovered in its disk (|Donati et all 
2005). There are no attempts known to date, however, 
to measure the magnetic field of EX Lupi. The inner rim 
of the high accretion region we measured (0.04-0.05 AU) 
translates to a period of 4-5 days assuming Keplerian 
rotation. The typical rotatio n time for T Tauri stars 
is a f ew days to a week (e.g. iRodriguez-Ledesma et al.1 
2009). The rotation period of EX Lupi itself is however 
not known yet. 

A self-regulated thermal instability theory assumes rel- 
atively high, constant inflow of mass (1O _6 M0 yr" 1 - 
1O" 7 M yr" 1 , IClarke et alJl989tlBeU fc TMT9M) . The 
constant mass transportation is hindered at a certain ra- 
dius, where the viscosity is simply too low to drain the 
incoming mass in time. The material piles up until the 
surface density becomes high enough to thermally ionize 
the gas. Near the ionization front, the disk opacity in- 
creases with temperature d ue to negative hydrogen ions 
addin g continuum opacity (|Faulkner. Lin, fc P apaloizou 
1983). Thermal instability is triggered by the runaway 
local heating of the disk. The transition front propagates 
outward until the instability is suppressed at the radius 
where the viscous heating is no longer effective to keep 
the disk warm. The radius at which the thermal insta- 
bility is turned off is a function of the mass accretion. 



It is 0.1 AU for a fiducial star of 1 M Q with 3 Rc^ with 
a con stant mass accretion of 3 x 10 _6 Mq (jBell fc LinI 
[1991 . 

Self-regulating thermal instability is well studied in the 
parameter space suited to FUors. The mechanism seems 
to work with accretion rates from M = 10~ 7 M Q yr -1 
to M = 10" 5 M Q yr -1 , and is able to reproduce 
the observed flaring magnitude at optical wavelengths 
(p»5 mag), the rise time to the maximum (~ 1 year), 
and the duration of the outburst (decades to a century) . 
These physical parameters are one order of magnitude 
smaller in the case of EX Lupi, with the quiescent phase 
accretion rate being M = 4 x 10~ 10 M Q yr -1 (|Sipos et al.1 

I2009f) which goes up to M = 2xlO~ 7 Af yr" 1 in the out- 
burst (jJuhasz et al.1120 j~0h with the rise time of a month, 
and duration of about a year. It is yet to be seen if the 
self-regulating thermal instability works for the small ac- 
cretion rate with no external triggers to cooperate. 

An interesting possible triggering mechanism is an em- 
bedded planet. When a protoplanet opens a gap, the 
mass transportation stops at its outer edge. The mass 
is banked up until the thermal instability is triggered 
(jLodato fc Clarkdl20?)l IClarke fc Sverl[l99l . With the 
additional degrees of freedom of the mass of the planet 
and its location, this model could cover a wider parame- 
ter space, and closely reproduce the observed properties 
of FUor outbursts. The triggering via embedded planet 
makes the instability propagate initially outside in, and 
qualitatively reduces the rise time. The smaller the mass 
of the planet, the smaller the outburst, triggered by the 
smaller accretion rate. It is still to be seen though, if the 
mechanism works for EX Lupi, as the exact parameter 
sp ace of EX Lupi outburst is not covered in the models 
bv lLodato fc Clarke! (pOfil . 

4.4. Variation of Line Profile Asymmetry 

Nevertheless, a possible planet in the disk is interesting 
in connection with the temporal variation of the line sym- 
metry in CO emission spectra (Fig. [5]). The line asym- 
metry indicates that the receding side of the disk was 
brighter on August 7, but by August 14 the approaching 
side became brighter. The receding part is again brighter 
3 days later, which again moved to the other side in 4 
days. The modulation of the line profile on the timescale 
of 3-4 days was k nown in the visible spectra of F U Or i 
and V1057 Cyg. iHerbig. Petrov. fc Duemnflerl (f200l 
discussed two possible origins of the variation, that the 
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modulation comes from the stellar spot and that it comes 
from the hot spot in the circumstellar disk, which may 
or may not be related to the hot accretion columns rising 
out of the disk. In the case of V1057 Cyg, the line mod- 
ulation is not accompanied by photospheric variability, 
which makes the stellar spot hypothesis difficult. On the 
other hand, the line modulation in V1057 Cyg is stable 
over 3 yrs, which also makes the hot spot on a differ- 
entially rotating disk unlikely. The question of whether 
the hot spot resides on the photosphere or on the disk is 
therefore still open. 

The period of the spectroscopic modulation of EX Lupi 
is hard to constrain with only six epochs of observations. 
If we simply take 7 and 3.5 days as half the period of 
the orbit in the first and the second half of August, the 
locations of the hot spot are 0.06 and 0.04 AU, respec- 
tively. This is within the span of the active disk that we 
derived from the outburst model, with the latter loca- 
tion close to the inner edge. Further spectroscopic mon- 
itoring with finer temporal sampling should be done in 
order to tell whether this represents real migration of a 
hot spot in the disk, as well as the theoretical study on 
the fate of protoplantear y clumps that fall into the star 
(jVorobyov fc Basul 120061 ). A more general discussion on 
the line asymm etry caused by a pla net embedded in the 
disk is found in lRegalv et all (|2010l ). 

5. SUMMARY 

Here is the summary of our results, and possible direc- 
tions for future investigations. 

• We reported the results of a spectroscopic moni- 
toring campaign of EX Lupi in the CO vibrational 
fundamental band at 4.6-5.0 /im. The observations 
covered six epochs in the outburst, which started 
in early 2008. 

• We detected high temperature CO gas in the 
emission-line spectra, that apparently consist of 
two distinct components. The equivalent width of 
the broad line emission (FWZI~ 150 km s -1 ) de- 
cays with the visual brightness of the star. The gas 
is vibrationally hot with the highly excited lines de- 
tected up to v=6. 

• The line profiles of the broad line component are 
double peaked, indicating that the gas is in a cir- 
cumstellar disk. We used a simple slab model to lo- 
cate the hot gas to the annulus 0.04-0.4 AU. There 
is strong asymmetry in the line profiles that varies 
with time on the order of days to week, implying 
the presence of a hot spot on the disk. 

• There is a second component of the line emission 
near the systemic velocity of EX Lupi with much 
narrower line profiles (FWHM«50 km s _1 ). The 
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narrow line emission is vibrationally cool; it exclu- 
sively consists of v=l-0 transitions. The line width 
indicates that the characteristic radius of the emit- 
ting gas is 0.4 AU, larger than the radius where the 
broad line emission arises. The equivalent widths 
of the narrow lines are constant over the whole ob- 
serving period. 

• An additional disk wind is found as absorption lines 
centered at —80 km s^ 1 toward the end of the out- 
burst. A similar disk wind was observed previously 
in V1647 Ori. 

• The two emission components are distinct in the 
decay rate, the excitation state, and the kinemat- 
ics. If the broad line emission traces the high ac- 
cretion region, its compact size favors local disk 
instabilities where a disk goes through bimodal ac- 
cretion. Major local disk instability theories, such 
as thermal instability, magnetohydrodynamical in- 
stability, magnetospheric truncation of a disk par- 
tially fit to the profile of the outburst, albeit our 
data did not allow to distinguish between these the- 
ories. The origin of the outburst involving global 
disk instabilities are unlikely. 

Most of the disk instability theories today have been 
tested against FUor outbursts, but little theoretical work 
has been done for EXors. It is not clear at the moment 
if the missing proper model for EX Lupi outbursts is due 
to an intrinsic difficulty in reproducing such small and 
frequent outbursts with low accretion rate, or is due to 
a lack of exploration in the parameter space that fits to 
EXor outbursts. Our attempts to remedy this are also 
hindered by missing fundamental parameters of EX Lupi. 
The period of the stellar rotation and magnetism, as well 
as an accurate simultaneous photometric calibration of 
the spectra at the time of the observation would help to 
improve the modeling of the outburst. 
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